addpath('npy-matlab-master\npy-matlab')
A1 = readNPY('map_LD_633nm_30K_NiI2_2_y.npy')
surf(A1,'edgecolor','none')
set(gca,'xlim',[1 21],'ylim',[1 21])
view(2)
%colormap(winter)
caxis([3, 12])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
%colorbar()
set(gcf,'Position',[0, 0, 300, 300])